Soliton model for proton conductivity in Langmuir films 
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5_( ' A soliton model for proton conductivity in Langmuir films is presented. The model 

Oh, contains three real scalar fields describing the hydrogen involved in the conduction, the hy- 

■ drophilic head of the Langmuir film, and the water. Soliton solutions that describe proton 
motion along the hydrogen bonds are found. Under compression of the film, the distance 
between the minima of the proton potential and the strength of the hydrogen bonds be- 
tween the film molecule and the water are changed. Such changes increase the probability 
of soliton creation. The model presented allows proton conductivity data in Langmuir films 

. to be explained. 

Or 

S 

O 

C/3 ■ 

^ , Research on the characterization of Langmuir films has led to the identification of a critical packing density, 

■ below which some properties of the system change significantly In particular, measurements of lateral 
conductance in Langmuir films suggested the authors of Ref . j|, " ] to postulate the conductance arises from 
proton transfer via the presence of hydrogen bonds in the film network. The likelihood of this assumption 

Oh is supported by the fact that in most of such films there is no room for ionic or electronic conduction. 

There are few theoretical efforts to understand the mechanisms behind such protonic conductance; see 
Ref. and references therein. In Ref. [|lO| the issue is addressed via a simple unidimensional mechanism 
^ , for proton transfer in the hydrogen bonds. In this model, the two oxygen atoms involved in the proton 
transfer are treated equivalently. This is an approximation, because one oxygen belongs to a water molecule, 
\ and the other to a hydrophilic headgroup. 

■ In the present work another model is proposed based on the soliton models of Refs. (ll- 13|. The model of 

o 



CZ3 



Ref. [|11| deals with proton conduction in quasi-unidimensional networks that arise due to hydrogen bonds. 
This model corresponds to a two-component chain, the proton subchain and the heavy nuclei subchain 
connected by the hydrogen bonds. It describes two real scalar fields, = (f)(x,t) and x = x(2^j^)i iii two- 
dimensional space-time. These fields correspond to the lighter {(p) and heavier (%) subchains, respectively. 
Q ■ The interactions between the two sublattices is described by the potential U — U{(t>, x), a functional of the 

■ fields. 

A related problem to the two-component chain is the diatomic one, which is modeled in Ref. ||l^. The 
' model introduces a derivative coupling between the two elements that form the diatomic chain. This model 
was further investigated in Ref. [^, where ideas of are combined to incorporate the derivative 

^ coupling between the lighter and heavier subchains of the two-component model. It was used to investigate 
• ^ , hydrogen-bonded systems. The model describes quasi-unidimensional systems that can be represented by 

■ repetitions of the basic entity X — H ■ ■ ■, where H and X represent the proton and the heavy nuclei, and 
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— and ■ • • the covalent and the hydrogen bond, respectively. The simplest example is water |T^, where X 
stands for OH . We can think of one-component and two-component models, but the two-component model 
is more realistic, because the network is visualized as a diatomic chain, composed of two sublattices, the 
lighter sublattice describing the protons, and the heavier sublattice the heavier component of the system. 
Langmuir films are spread and compressed on an aqueous subphase (substrate in our case), and we recall 
that the experimental measurements in such films are in general done at room temperature. The film is 
considered an isolated system, having no interference from the water substrate. The water substrate only 
provides hydration water that participates in the H-bond network. In the model, we do not have to assusme 
any specific crystaline structure for the hydrated water. 

The main goal of the present work is to offer a new model to understand the proton conductivity problem 
in Langmuir films. This model appears through the hydrogen bonds that connect the hydrophilic head of the 
film molecules to the water, leading to a mechanism that allows the presence of solitons, that drive the proton 



motion in the network. To do this, we first extend the model of Ref. |13| to the case of a three-component 
chain, which seems to be more appropriate to the investigation of such films. The three-component network 
in the Langmuir films can be schematically represented by repetitions oiX~H---Y — H---. It is composed 
of protons and the X and Y groups that represent OH and the amphiphilic molecule that characterizes the 
film, respectively. They are spatially represented in Fig. 1 in the case of the aliphatic acid. 

We justify a three-component model for such system recalling that even when one neglects the tail con- 
tribution to the motion of the head, it is much heavier than the water counterpart. The proton subchain is 
much lighter than the other two subchains, so we consider the lighter sublattice coupled with each one of 
the two other sublattices. We assume the two heavier sublattices do not interact with one another by any 
other mechanism. That is, the film molecules only interaction with the water is via the hydrogen bonds. 

Our model uses three real scalar fields 0, x, and t/", in bidimensional space-time. 4> = 4>{x,t) describes 
translational motion of protons, and x = xi^i ^) ^^'^ "0 = "0(2^1 1) the motion of OH groups and amphiphilic 
molecule subchains, respectively. We use standard notation, with a;" — {t,x), Xa — {t, — x), and h = c = 1. 
The Lagrangian density is given by 

= ^ da(j)d°'<j) + ^ Ail daXd^X + ^ M2 d^'ipd^'tp 

The parameters /ii, /12, and vi, 1^2 are real and positive, and U — U{<f>^ Xi ■0) is the potential. Here we consider 
the protons as unit mass particles, so the parameters /zi and /i2 can represent the mass ratio between the 
OH group and the hydrogen, and the amphiphilic molecule and the hydrogen. The other two parameters 
vi and V2 describe the derivative coupling between the proton and water and the proton and amphiphilic 
molecule, respectively. 



The equations of motion are 

dx'^ ' ^ dx'^ ' ^ dx^ 



d^<i> , d^x ^ d^i> dU 

'^i 3-2 +^23-7 = -TTT, (2) 



d^X^ dU 

and 

^ ^ <l± _d£_ 

dx'^ dx^ 

We consider the case in which the potential is defined by some smooth function W = W{(t), x, in the form 

u{4>.x,i')^\wl + \wl^]^wl (5) 
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stands for dW/d(j), and so forth. This specific form of the potential, together with ideas first introduced 
in Refs. [p"5|-|l7| have been recently used to model topological twistons in crystalline polyethylene 18 19[|. 
Other investigations related to the subject can be found in Refs. 120-12^ and in references therein. 




FIG. 1. Schematic view of the three-component network, (a) The proton, the OH group, and the head of the film, 
are represented by small, medium and large circles, respectively, (b) Shows explicitly the atoms of the diagram (a). 
The hydrophobic tail of the amphiphilic molecule is represented by the gray area. 

In Ref. the derivative coupling between the two heavier subchains was considered in the two-component 
model there studied. In the present work we are considering the three-component model just introduced, 
and no coupling between the two heavier subchains is taken into account, since they seem to couple mainly 
through the hydrogen bonds. We follow this reasoning and neglect further couplings between the x ^^'^ V' 
fields in the potential. Also the complexity of the problem can be reduced by assuming the heavier subchains 
contribute insignificantly to the potential, at least when compared to the lighter subchain. Therefore, in this 
approximation U = U{(j)) depends only on the mobility of the proton field, which yields to 

U{'P) = \wl. (6) 
The equations of motion for static field configurations are reduced to 

(P(h (Py (Pi' , s 
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and 



^'^ dx^ ^ dx^ 



0. 



(8) 



(9) 



We consider fields whose derivatives present similar asymptotic behavior. In this case the above equations 
(0) and (0) change to, after setting to zero the integration constants, 



dx 
dx 
dip_ 
dx 



vi d(j) 
dx ' 
1^2 d(j) 
^2 dx' 



(10) 
(11) 



For static configurations the non- vanishing components of the energy-momentum tensor T°'^ can be written 
as T^^ = Tg + U{(t)) and T^^ =Tg — U {(/)), where Tg represents the gradient contribution. This is given by 



d<j) 



dx 



1- +Mi 3- 



dx 



dx 



-2vi 



d4>\ ( dx 



dx I \dx 



2V2 



dip 



d(j}\ f dip 



dx I \ dx 



We use the solutions ( p^ ) and (|T^) to write 



rpW _ ^ 

2 

rpll _ }_ 



B 



2 f^\\f^ 



\dx ) ^ \ 



d(P 
dx 



/dW 
\~d^ 



where B^ = 1 — vf/ni — v\j y^2- B is considered real and positive. 

The energy density of the static solutions e{x) is identified with T^", and it can be written as 



(12) 



(13) 
(14) 



e^-[ B---W, 



dx 



To minimize the energy, the first term in the above expression is set to zero 

dW 



d(j) 
dx 



with minimum energy solutions given by 



E ^B\/\Wl 



(15) 



(16) 



(17) 



where ISW = Vl^[0(cx))] — M^[(/)(— cxd)]. Solutions to the first-order equation ( p_6[ ) are known as Bogomol'nyi- 
Prasad-Sommerfeld or BPS solutions ||2^,^. We notice that T^^ vanishes for the BPS solutions, as expected. 
For the proton self-interaction we consider 



W{cP) = iA03 _ ^^2^ 



(18) 
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which gives the desirable double-well potential 

U{cP) = ^X'{<j>'-A'r, (19) 

where A is a real and positive dimensionless parameter, A is also real and has dimension length"^. has 
minima located at (j)± = ± A and the barrier height A'^/2. Since the parameter A gives the minima of the 
potential, it is directly related to the distance between neighbor oxygens in the film network. This allows 
the presence of two degenerate ground states, one with all the protons at the position A, and the other with 
the protons at —A. These two states can be represented by repetitions of the basic entities 

X-H • •• Y- H • •• 

X ■■■ H-Y ■■■ H- 

The presence of two degenerate minima allows the appearance of solitons, which are extended solutions with 
finite energy that connect the two degenerate ground states. 

To investigate the presence of soliton solutions, we notice that in the present model the first-order equation 
( p^ ) becomes 

^ = XB-\c^^ - A^), (20) 
dx 

It is solved to give the kink solutions 

(t>{x) ^ A tanh {AB-^ x) . (21) 

Here x stands for Xx, and is dimensionless. The parameters A and B determine the energy {A/3)A^^ B^^ , 
and the width I ^ B/A oi the soliton solutions. A is related to the distance between neighbor oxygens in 
the H-bond network. It is an important parameter because under compression, conductivity and surface 
potential of the film change significantly when the critical area per head group limit is reached [^J-Q. B is 
an effective parameter, that depends on the fundamental parameters fii, fj,2,vi, V2 in a specific way, imposed 
by the field-theoretical model here considered. Although our model is microscopic, we can think of B as a 
phenomenological parameter, used to infer the width and the energy of the soliton solutions. The soliton 
solutions spring in responce to nonlinear effects, and describe proton mobility in the network. The above 
solutions show that protons migrate from A to —A, but this migration takes several units in the chain 
network, given in accordance with the width / B/A of the solution. The main characteristics of the 
solution is illustrated in Fig. 2, where the soliton and the asymptotic values in the chain network is depicted. 
Evidently, under the presence of an external electric field the soliton may move, giving rise to a steady current 
along the direction dictated by the external electric field. The picture is similar to that in the original work 
on solitons in hydrogen-bonded network [0 . 

We consider a spatial arrangement of the molecules as in Ref. As the film is compressed, the area 

per amphiphilic unit reaches a critical limit, in which the distance between two consecutive COOH groups 
is around 7 A; see Fig. 1. At this point, the oxygens from the OH and film head groups are about 2.4 
A apart. The hydrogen bonds get stronger, bridging the water to the film head tightly, which forms the 
network that allows the proton conduction. Experimental data shows that the conductivity increases with 
decreasing of film area. This aspect can be incorporated in the model if we assume the energy of the kink 
decreases with area per amphiphilic molecule. We do this by requiring the parameter A to decrease with 
the distance between the head groups. 
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FIG. 2. Schematic view of the sohton solution, showing the corresponding amphtude A and width I ~ B/A. We 
also depict the asymptotic spots in the chain network, to illustrate that the soliton is localized in a finite region, 
involving several units of the hydrogen-bonded chain. 

When A decreases, the parameters i>i and 1/2 that control the derivative couphng should increase, because 
when the head of the amphiphilic molecule gets closer to the water molecule, the derivative coupling should 
become more effective. Since B is given by i? = 1 — vf/ fxi — i'2/fJ-2, for increasing Vi and 1/2, it should 
consequently decrease. Although we do not know the exact way B decreases, we can compare it with A. 
We can consider, for instance, the simplest possibility in which the ratio A/B remains constant. This 
approximation introduces two important consequences: i) the soliton width depends on A/B, so it should 
not depend on the film compression; ii) the energy of the soliton depends on A^/B, so it should only 
depend on A^. With this assumption, our model predicts that the energy of the soliton should vary at 
the rate E/Eq = A^/Aq, for some reference value Aq. As the energy is reduced, the probability of soliton 
creation is increased, yielding a higher proton conductivity in agreement with the experimental results. In 
summary, the formalism presented here allows proton conductivity data in Langmuir films to be explained. 
Because the physical grounds have been established, one can now extend the soliton model to treat the data 
quantitatively, which will be the subject of our further investigation. In addition, the framework may be 
applied to proton conductance in more involved systems such as cell membranes, where proton conductance 
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is believed to play an important role Also it may be useful in the attempts to investigate effects from 
alcohols on lateral conductance ||2^ ]. 
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